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Quasi-thermal Universe and its implications for gravitino produciton, baryogenesis 

and dark matter 
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Under general circumstances full thermal equilibrium may not be established for a long period 
after perturbative or non-perturbative decay of the inflaton has completed. One can instead have a 
distribution of particles which is in kinetic equilibrium and evolves adiabatically during this period. 
Number- violating reactions which are required to establish chemical equilibrium can become efficient 
only at much later times. We highlight some of the striking consequences of such a quasi-thermal 
Universe. In particular, thermal gravitino production yields no bound on the maximum temperature 
of the primordial thermal bath alone. As another consequence, the lower bound on the mass of the 
lightest right-handed (s)neutrino from thermal leptogenesis can be ^> 10 9 GeV. Depending on the 
phase, a Wino or Higgsino considerably lighter than TeV, or a Bino in the bulk region, can be a 
viable thermal dark matter candidate. Finally the electroweak symmetry may never be restored in 
the early Universe, therefore weakening any hopes of realizing a successful electroweak baryogenesis. 



The detection of the cosmic microwave background 
anisotropies and the fluctuations on scales larger than 
the size of the Hubble radius at the recombination era Jl[ 
is a strong indication of an early inflationary epoch [2j. 
An important aspect of inflation is the reheating, which 
starts with the inflaton decay and, eventually results in a 
thermal bath of elementary particles. The energy density, 
p, and the number density, n, of a distribution of rela- 
tivistic particles at a full thermal equilibrium are given 
by 



p = g* (vr 2 /30) T 4 , n = (C(3)/tt 2 ) T 3 



(1) 



where T is the temperature and g* is the number of de- 
grees of freedom. The average energy of particles is given 
by: (E) = (p/n) ~ 3T ~ p 1 / 4 . 

The inflaton decay can be a perturbative [3| , or a non- 
perturbative Q process, while in the models of running 
mass inflation the inflaton condensate fragments into un- 
stable lumps slowly decaying through surface evapora- 
tion In all these cases the inflaton decay results in 
a distribution of particles which is far from full thermal 
equilibrium. 

Let us denote the inflaton mass (close to the minimum 
of its potential) and the decay rate by rn^ and Td respec- 
tively. At the end of the inflaton decay: p « 3 (Fd-Mp) 2 , 
where Mp = 2.4 x 10 18 GeV is the reduced Planck mass. 
When the inflaton decays perturbatively, Td -C m 2 /Mp 

and (E) w which results in (E) ^> p 1 / 4 . On 

the other hand for the non-perturbative inflaton de- 
cay as in the case of preheating, Td ~ 10 and 
(E) > 1Oto0 0], implying that (E) < p 1/4 . In both 
the cases reaching full equilibrium requires that (E) and, 
from the conservation of energy, n change and settle to 
their respective equilibrium values given by Eq. |T]). 

Reactions which conserve the total number of particles 
re-distribute the energy and lead to kinetic equilibrium, 
as well as relative chemical equilibrium, among differ- 
ent degrees of freedom. Full thermal equilibrium also re- 
quires number-violating reactions, which are higher order 



processes, so that the chemical equilibrium is achieved 
for all species. An important point is that in general the 
inflaton decay, number-conserving reactions and number- 
violating reactions occur at (very) different rates. In 
consequence, kinetic and chemical equilibrium can be es- 
tablished at different time scales long after the inflaton 
decays. 

Indeed this is the case after preheating: the inflaton de- 
cay products (and the fields to which they are coupled) 
reach kinetic equilibrium very rapidly, but full equilib- 
rium is established over a much longer period jfj . In the 
case of perturbative inflaton decay, the 2 — ► 2 and 2^3 
scatterings with gauge boson exchange have the main 
role in building kinetic and chemical equilibrium respec- 
tively Q. These scatterings are therefore suppressed if 
the relevant gauge bosons have large masses in the early 
Universe. 

This generically happens in supersymmetric theories. 
Supersymmetry guarantees flat directions (made up of 
squark, slepton and Higgs fields) which acquire a large 
vacuum expectation value (VEV) during inflation |8j. 
Flat directions with a large VEV usually enhance the 
inflaton decay [9]. However, similar to the Higgs mech- 
anism, the flat direction VEV (denoted by (ip)) sponta- 
neously breaks gauge symmetries and induces a mass m g 
to the gauge bosons, where m 2 ~ a(ip) 2 (a is a gauge 
fine structure constant). Note that m g has no bearings 
at low energies/late times since it vanishes once the flat 
direction settles to its origin 1 . 



1 Particles which are coupled to the flat direction acquire a large 
mass oc (</?}, and hence promptly decay to massless particles (i.e. 
the ones which are not coupled to the flat direction). Such a 
decay can increase the number of particles by a factor of 2 or 3, 
depending on whether the dominant channel is 1 — > 2 or 1 — > 3 
decay. However this alone is too small to lead to full chemical 
equilibrium, as it needs the number of particles to increase by a 
huge factor (typically > 10 3 ) after perturbative inflaton decay. 
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Kinetic equilibrium is usually built at a rate I\i n ~ 
a 2 n/m 2 , while full equilibrium is established at a rate 
I\hr ~ a 3 n/i7ig (n is the number density of infla- 
ton decay products) 0]. For sufficiently large (ip) ~ 
O(10 13 - 10 18 ) GeV and a ~ 1CT 2 , it turns out that 
Tthr *C Td- This substantially delays thermalization and 
results in a (very) low reheat temperatures (for example, 
T R ~ O(TeV) if (tp) ~ 10 18 GeV) 0. 

Therefore under general circumstances the Universe 
enters a quasi-thermal phase. During this phase the co- 
moving number density and comoving average energy of 
particles remain constant. It is then paramount that 
the Universe reaches full equilibrium before nucleosyn- 
thesis i.e., Tr ^ O(MeV), which is known to toler- 
ate hardly any departure from thermal equilibrium. The 
aim of this paper is to show striking consequences of such 
a quasi-thermal epoch, such as particle production and 
phase transitions. 

During the quasi-thermal phase we can parameterize 
the energy density and the number density of the i— th 
degree of freedom as: 



Pl = C t A t (tt 2 /30) T 4 , n t = Ai (C(3)/tt 2 ) T 3 



(2) 



where Ci, Ai ^ 1 indicate deviation from full equilibrium 
(i.e., Ci = Ai = 1 for bosons and A4 = 3/4, C* = 7/6 
for fermions). Note that relative chemical equilibrium 
among different species may be established much later 
than kinetic equilibrium (as it happens after the pertur- 
bative inflaton decay [l(|), and hence have different 
values in general. If Ai < 1, the number of particles is 
smaller than that in full thermal equilibrium. We can 
then assign a negative chemical potential /ij to the i— th 
particle where Ai — exp(fii/T), and Ci ~ (E)/3T as 1. 
This is a typical situation for the perturbative inflaton 
decay 2 . 

The situation will be more complicated for Ai > 1, be- 
cause the number of particles is (much) larger than the 
full equilibrium value, thus implying a positive chemical 
potential. For bosons with a typical mass m, < T, we 
have Hi < m^. This implies population of the ground 
state and renders large occupation numbers. Once n,; 3> 
T A /rrii, the ground state dominates both the number 
and energy densities. We can then use Eq. ([2]) where 
Ai ^> 1 and Ci ~ rrii/T <C 1. Similarly, for a degener- 
ate Fermi gas with fa 3> T, Eq. is applicable with 
Ai ~ (fa/T) 3 > 1 and C l ~ fa/T > 1. At the end 
of preheating the average kinetic energy and the (effec- 
tive) mass of particles in the plasma are approximately 



2 Since in this case particles are produced in one-particle decay of 
the inflaton, the distribution is peaked around the average en- 
ergy. The 2 — > 2 scatterings (which become efficient shortly be- 
fore complete thermalization) smooth out the distribution. This 
implies that the reheat plasma is close to kinetic equilibrium but 
grossly deviated from chemical equilibrium, hence is in a quasi- 
thermal state. 



equal This implies that ~ T and Cj ~ 1. There- 
fore, in cases of physical interest, we can take Ci ~ 1 in 
Eq. ©. 

The quasi-thermal phase then lasts until H ~ Tthr, 
where r t hr is the rate for number-violating reactions. 
During this period Ai remain practically constant and 
T cx i? -1 , rii oc R~ 3 , pi oc R~ 4 (with R being the scale 
factor of the Universe). This implies an adiabatic evolu- 
tion of the plasma during which 



H ~ A 1/2 (T 2 /3M P ) , A = ^2Ai 



(3) 



At full equilibrium A = g*, where ~ 200 in the 
minimal supersymmetric standard model (MSSM). The 
plasma has its maximum temperature T max at the end of 
the inflaton decay, i.e., when H ~ T^. From Eq. we 
obtain, 



A ~ 10 



r^Mp/T, 



1. 1 : ax ) 



(4) 



Note that fa and Ai are determined by the branching ra- 
tio for the inflaton decay to the i— th degree of freedom. 
Perturbative inflaton decay generically yields Ai <C 1 
(thus A <c <?*)■ For example, consider a gravitation- 
ally decaying inflaton Td ~ rrA/Mp with = 10 13 
GeV. This leads to T max - 10 12 GeV (note that (E) ~ 
3T max as m<f,/2) and A ~ 10~ 8 . On the other hand, pre- 
heating results in the opposite situation where Ai 1 
(and i > j,). As mentioned earlier T max ~ lOm^ and 
T d ~ lO~ 2 m in this case 0. Then, for m = 10 13 GeV, 
this yields A ~ 10 4 . 

Once H ~ Tthr, at which time the temperature has 
dropped to Tp, the number- violating reactions become 
important and lead to full equilibrium (i.e., Ai ~ 1) very 
quickly. 

We define the reheat temperature Tr as the tem- 
perature of the distribution after the establishment of 

1/2 

full thermal equilibrium, therefore Tr ~ (Tthr-Mp) 
In general we can have r t hr <S Td, and hence Tr <C 
(Fd-Mp) 1 / 2 , which is often used in the literature Q. The 
conservation of energy implies that Tr = (A/g*) 1 ^ 4 Tp. 
The entropy density, which should be used to normalize 
a given number density, is given by 



= A 3 ' 4 gl /4 (2tt 2 /45) T, 



(5) 



Let us now consider the production of a (stable or long- 
lived) particle x with mass m x in the quasi-thermal 
phase. If the couplings of x to the fields in the thermal 
bath are sufficiently small, it will never reach equilibrium. 
Then its abundance will be subject to only production 
and redshift, and we will have 



3Hn, 



E 



A,A, 



(C(3)A 2 )>> cl ) T e 



(6) 



Here denotes the cross-section for production of x 
from scatterings of the i— th and j— th particles and the 
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sum is over all particles which participate in x produc- 
tion. Also () and the dot denote thermal averaging and 
proper time derivative, respectively, and we have used 
Eq. |J2J|. Note that \ production is kinematically possi- 
ble in the interval m x <J T < T max . 

One comment is in order before proceeding further. 
In supersymmetric theories the last stage of inflaton de- 
cay which reheats the Universe is generically perturba- 
tive and results in Ai < 1 [10]. In this (physically rel- 
evant) case Eq. ^ can be relibaly used to calculate % 
production in the quasi-thermal phase. However the sit- 
uation will be more subtle when A; 3> 1. In this case 
we have a dense plasma where the occupation number 
of particles is large (i.e. ^> 1 for bosons and ~ 1 for 
fermions). Such large occupation numbers, besides fi- 
nal state effects which enhance or suppress the interac- 
tion rates, result in large induced masses and significant 
off-shell effects. Therefore the situation is considerably 
more complicated and a pro per (non-equilibrium) quan- 
tum field theory treatment [12j | will be required in this 
case. Thus the following expressions of particle abun- 
dances must be taken with care and considered only as 
order of magnitude estimates when Ai 3> 1. 

Gravitino- If \ is the gravitino, ~ a/Mp. In gen- 
eral, an exact calculation of the gravitino abundance is 
complicated since various scattering processes have dif- 
ferent cross-sections (similar argument will be followed 
for the rest of the paper). The result now depends 
on the values of Ai, and hence is model-dependent. 
The maximum abundance is obtained when the plasma 
mainly consists of gluons as gluon + gluon — > gluino + 
gravitino scattering has the largest cross-section. Fol- 
lowing the standard calculations [13] ]. and with the help 
of Eqs. (|3l5p . we then find an upper bound (for a more 
detailed calculation see [lOj]) 



(n 3/2 /s) ~ A 3 / 4 (T max /10 10 GeV) x HT 12 



(J) 



in the case of MSSM. For a TeV scale gravitino the suc- 
cess of BBN requires that 14 1 



T max < A" 3 / 4 x 10 10 GeV . 



(8) 



This bound depends crucially on the parameter A. To 
illustrate, let us again consider the example of a gravi- 
tationally decaying inflaton for which A ~ 10~ 8 . It is 
interesting that despite having T max ~ 10 12 GeV, the 
constraint in Eq. ||5J) is comfortably satisfied. Note that 
now the bounds from thermal gravitino production on 
the inflaton sector can be substantially relaxed. There- 
fore within the regime of perturbative decay there is no 
strict bound on maximum temperature of the Universe. 

On the other hand in the case of preheating (where 
A ^S> 1), the constraint from the overproduction of grav- 
itinos become: < T max <C 10 10 GeV. For such a small 
inflaton mass chaotic type inflation will be a disaster from 
the point of view of matching the observed temperature 
anisotropy, one would have to resort for alternative sce- 
narios as discussed in Ref. (l5j . Otherwise one would 



require a larger gravitino mass which would decay early 
enough as discussed in Ref. [l6[ . In passing we note that 
gravitinos are also created during preheating from vac- 
uum fluctuations [l7| . This non-perturbative production 
has nothing to do with thermal scatterings, and hence is 
not modified in the presence of a quasi-thermal phase. 

Right-handed (s)neutrino- If \ is the lightest right- 
handed (RH) (s)neutrino, its decay can generate the ob- 
served baryon asymmetry of the Universe. In the thermal 
leptogenesis scenario these (s)neutrinos are mainly pro- 
duced from scatterings of the left-handed (LH) (s)leptons 
and the top (s)quarks [l8|. The exact result is again 
model-dependent. However, the generated baryon asym- 
metry in the quasi-thermal phase is maximum when the 
plasma mainly consists of the LH (s)leptons and the top 
(s)quarks. Applying Eqs. (|3I5[) to the standard calcula- 
tions ,18[, then results in an upper bound (for a more 
detailed calculation see (Toj'): 



r/ B = n B /s ~ A 3/4 k (Mi/10 9 GeV) x 3 • 10" 



(9) 



where Mi is the mass of the lightest RH (s)neutrino, 
and k < 10 _1 is the efficiency factor accounting for the 
processes involving the RH states. Generating sufficient 
asymmetry 77B — 0.9 x 10~ 10 then sets a lower bound on 
A 3 / 4 Mi instead of Mi alone: 



Mi > A~ 3/i x 10 9 GeV . 



(10) 



For a perturbative inflaton decay (^4 <C 1), this results 
in Mi > 10 9 GeV. However, Mi <, 3T max w m /2 is 
also required in order to kinematically produce the RH 
(s)neutrinos, where < 10 13 GeV. This, together with 
Eq. (p~0|) . leads to a lower bound A > 10~ 5 . This implies 
that for a gravitationally decaying inflaton A ~ 10~ 8 is 
too small to allow successful leptogenesis. 

In the case of preheating, where A 3> 1, Eq. (fit)]) then 
implies that successful leptogenesis is possible for much 
lighter RH (s) neutrinos Mi C 10 9 GeV. One point to 
note is that the gravitino abundance has the same de- 
pendence on A as the baryon asymmetry, see Eq. (|7l9p . 
Therefore thermal leptogenesis remains marginally com- 
patible with thermal gravitino production in a quasi- 
thermal Universe as well. 

Thermal dark matter- If x is one of the particles in the 
thermal bath, its distribution will follow Eq. (|2|). This 
happens, for example, when x is t ne lightest supersym- 
metric particle (LSP). If Tp < m x , the x quanta become 
non-relativistic before full thermal equilibrium is estab- 
lished. Again assuming that the plasma mainly consists 
of x, its equilibrium number density and energy density 
are in this case given by: 

n x ~ ^(m x T/27r) 3/2 exp(-m x /T) , p x = m x n x . (11) 

The annihilation of x quanta, happening at a rate r ann = 
(c x /m x ) n x , keeps them in equilibrium until the time 
This occurs at T = Tf, where the 
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freeze-out temperature, Tf , follows: 

m x /T { = x { , X{ ~ 0.51ilA + In [0.2 c x (M P /m x )] . (12) 

The entropy release upon full thermalization, given by 
Eq. |5]), results in the final abundance 

n x /s ~ A- 1/4 x t (0.2/c x ) (m x /Mp) . (13) 

We remind that in the case of full equilibrium, the uni- 
tarity bound implies that c x < 4n, and hence the ther- 
mal abundance exceeds the dark matter limit (n x /s) < 
5 x 1(T 10 (1 GeV/m x ) if m x > 100 TeV 0. Freeze- 
out during the quasi-thermal phase, see Eq. (fl~3|) . con- 
strains A _1 / 4 m x instead. Thermal dark matter with 
m x > 100 TeV will be now viable if A > 1. On the 
other hand, the upper bound on m x will be strengthened 
ifA<l. 

There are particularly interesting consequences for the 
LSP dark matter [2(|. In full equilibrium, viable Wino- 
or Higgsino-like dark matter requires that m x > 1 TeV. 
This is considered too heavy if weak scale supersymme- 
try is the solution to the hierarchy problem. Freeze-out 
in the quasi-thermal phase with A <C 1 allows m x to 
be (much) closer to the ALEPH bound on the lightest 
neutralino mass (54 GeV at 95% c.l. 21]). If full equilib- 
rium is assumed, a Bino-like LSP overdoses the Universe 
in most parts of the parameter space due to its much 
smaller annihilation cross-section. However, see Eq. (p~3|) . 
A ^> 1 can significantly extend the acceptable part of the 
parameter space to the bulk region as it lowers the LSP 
abundance. 

Electroweak phase transition- Finite temperature ef- 
fects make a typical contribution n/T to the (mass) 2 of 
every scalar field which is coupled to the thermal bath. In 
full equilibrium, Eq. |T]), thermal corrections are ~ +T 2 . 
Therefore, at sufficiently high temperatures, they can 
lead to restoration of a broken symmetry [22j ■ Thermal 
masses can also affect the dynamics of cosmological Bose 
condensates by triggering their early oscillations [231 ]. 

For a quasi-thermal distribution, Eq. ©, thermal cor- 
rections are ~ AT 2 and depending on A there will be dif- 
ferent situations. As a specific example, consider the elec- 
troweak symmetry breaking happening at a scale Wew ~ 



100 GeV. This implies that the electroweak phase tran- 
sition occurs at a temperature T cr ~ A _1 / 2 mEw- 

If A < (m EW /T max ) 2 , then T a > T max and the elec- 
troweak symmetry will be broken after inflation. If full 
thermalization occurs very late, such that < 100 GeV, 
it will never be restored and the electroweak baryogenesis 
will not be possible. 

On the other hand, T cr < 100 GeV if A > 1. Such a 
late phase transition will also be problematic for the elec- 
troweak baryogenesis, since sphaleron transitions are not 
active at that time. Note that A 3> 1 can also result in 
a large mass mw ^> T for the electroweak gauge bosons 
at T > T cr . The sphaleron interaction rate, which is 
oc exp (— mw/T), will then be exponentially suppressed 
even before the electroweak phase transition. 

Discussion- The most important message of this paper 
is that thermal particle production does not solely de- 
pend on the temperature of the Universe. The reason is 
that under general circumstances the early Universe can 
enter a quasi-thermal phase after the inflaton decay. This 
epoch, during which the Universe is far from full thermal 
equilibrium, can last very long and include a wide range 
of temperatures, particularly those which are relevant to 
the production of gravitinos, baryon asymmetry or dark 
matter. 

In particular, we showed that thermal gravitino pro- 
duction yields no bound on the maximum temperature 
of the Universe alone. The lower bound on the light- 
est RH (s)neutrino from thermal leptogenesis can also 
be very different from 10 9 GeV. Depending on the ther- 
mal history, a Wino- or Higgsino-like LSP much lighter 
than TeV, or a Bino-like LSP in the bulk region, can be 
a viable thermal dark matter candidate. A detailed ac- 
count of thermalization in supersymmetric models and 
its cosmological consequences will be given in a separate 
publication [lOj. 
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